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SUMMARY 


Variational methods for the approximate solution of subsonic and 
transonic flows of a compressible fluid before the occurrence of shock 
waves have been cai*ried out in previous papers. The methods fall as 
soon as the shock waves occur as the flow behind the shock waves now 
becomes rotational and has variable entropy. Since most transonic flows 
are accompanied by shoc]^ waves, a method which allows for shock waves 
and variable entropy is necessary for the study of such flows. By 
modifying Bateman’s variational principle for Irrotational flows, it is 
shown that a variational principle for flows with rotation and variable 
entropy can be obtained. By applying this variational principle to the 
regions of flow behind shock waves and Bateman’s original principle to 
the other regions in the flixld, shock equations can be directly obtained. 

A procedvure for computing numerical solutions for such flows is suggested, 
and a nximerlcal example is carried out. At hl^ Mach number above a 
certain limiting valiie, the restilts show that Irrotatiorial flow falls. 
However, by inserting shock waves and allowing a part of the flow to be 
rotational, computation indicates that solution exists again. 


INTRODUCTION 


In previous papers (references 1 to 7) the senior author and his 
associates have succeeded in applying the variational methods to the 
study of subsonic and transonic flows of a compressible fluid past 
arbitrary bodies before the occvtrrence of shock waves. Numerical examples 
were carried out in the case of the flow past a circular cylinder, an 
elliptical cylinder, a Kaplan bump, a sphere, and an ellipsoid. The 
results were found to check excellently with those computed by other 
methods. These results Indicate that the variational method will give 
good approximations to flows past either thick or thin bodies and at 
both low and hi^ Mach numbers. The method as formulated, however, can 
only be applied to irrotational flpws. As soon as shock waves occur the 
method fails because the flow behind the shocks then becomes rotational. 

As most transonic flows are accompanied by shock waves, different methods 
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whlcli allow the rotatlonality of the flow as well as variable entropy- 
must be formulated if any significant results are to be obtained from 
the studies of such flows. 

A reexamination of Bateman’s variational prin'. rples (reference 8) 
indicates that with some modification the principle in terms of the 
stream function which was originally formulated for irrotational 

flows, can also be applied to ro-bational flows with entropy change. The 
resulting -variational Integral is the same as the one recently obtained 
by Lin and Rubinov (reference 9)- In the study of transonic flows this 
integral is to be applied in the region of flow after shock waves and 
the original Bateman integral in the other regions. It can be sho-wn 
that the shock equations are directly obtainable from these principles. 

With the -variational principles ob-tained, a direct method for the approxi- 
ma-te solution of transonic flow with shock waves may again be formulated 
-following the Raylei^-Rltz procedure. The actml carrying out of such 
a method however was found too laborious. Instead, Galerkln’s method 
may be used which shortens the numerical work to a great extent. In 
this report approximate solution in -fche case of the transonic flow past 
a circular cylinder has been carried out. The refe’olts show that when 
the Mach num'ter increases -fco a cer-fcain limiting value without allowing 
for shock -waves the variational method does not ha-ve a solution. This 
probably Indlca-tes the bi-eakdo-wn of the irrotational flow. By allowing 
the occurrence of shock waves, solution again exists. 

This work was conducted at the Daniel Guggenheim School of Aeronautics, 
College of Engineering, New York University, under the sponsorship eind 
with the financial assis-fcance of the National Advisory Committee for 
Aeronautics. The authors are indebted to Professor F. K. Telchmann for 
his kind Interest. 


SYMBOLS 

a velocity of sound 

Ann^®mu:’^’mn^®*mn undetermined parameters 

Cp,cv specific heat at constant pressure and -\rolume, 

respectively 

CQ_,C 2 ,C 2 constants 

M Mach number (q/a) 


n 


normal distance 
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P 

q 

s 

u 

V 

7 

<p 

p 

ca 

U 

Subscripts 

x,y,r,0 

5 

O 


pressure 

magnitude of velocity vector, with components u 
and V 

entropy 

velocity component in x-direction 
velocity component in y-direction 
coordinates of a point in the fluid 
isentropic exponent |cp^ c^ 

velocity potential function 

stream function 

mass density 

rotation of flow 

velocity of undisturbed stream 

denote differentiation in corresponding direction 

isentropic stagnation conditions for undisturbed 
stream 

conditions of undisturbed stream 


GOVERNING DIFFERENTIAL EQUATIONS AND BOUNDARY CONDITIONS 


Consider a steady two-dimensional flow passing a cylinder of 
arbitrary shape. Most aerodynamic problems are concerned with either 
flow from rest or flow which is parallel and uniform at Infinity. In 
such cases before the occuirence of shock waves the flow is irrotational. 
As the speed of the flow increases, shock waves occur. For example, the 
flow past an airfoil is shown in figure 1. The flow is asstaned to be 
subsonic at infinity. Becaxise of the presence of the airfoil, if the 
Mach nimiber of the flow is stifflciently hi^ there will be a region of 
supersonic flow near the surface of the airfoil, as shown by the dotted 
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lines in the figure. Shock waves will occur at some points on the 
surface in these supersonic regions and will terminate in the flow 
where local Mach number is eqxial to 1 . Then in the region extending 
to infinity behind the shock waves and hounded hy the stream lines 
passing through the points where the shock terminates, the flow is 
rotational and with variable entropy. For convenience, this region 
will be denoted by D2, and the region outside D2, by D]^. Then 

in Dp the following equations must be satisfied: 

(1) The two equations of motion 

" "Px/P 

( 1 ) 


where u and v are the velocity con^jonents in the x- and y-directions, 
p the density, and p the pressure, and the subscripts x and y 
indicate partial differentiation in the corresponding direction. 

(2) The equation of continuity 


UVt, 


Wy = 


-Py/P 




J 


(pU)x + (pV)y = 0 


( 2 ) 


(3) The equation of state 


P = cpp 7 ( 3 ) 

where cp is a constant and 7 is the ratio of specific heats. If the 
subscript s is used to Indicate the stagnation conditions, cp = Ps^Ps^- 

(ll-) The irrotationality condition 




Uy - Vx = 0 
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With aid of equation (4), equation (l) can te integrated to give 


*(r - l)p 2 


<1 

+ ^ = c. 


(5) 


where q^ = u^ + and C2 is another constant. 

% 

The boundary conditions are that at infinity the velocity components 
u and V are equal to the given values^ and on the solid boundary of 
the airfoil the noimal component of the velocity vanishes. 

In region D2 the flow becomes now rotational and the entropy S 
is no longer a constant. Here the governing differential equations are 
as follows: Equation of motion (l) and the equation of continuity (2) 

are still valid. The equation of state has to be changed to the 
following fom: 


p = 036®/'^'^ p 7 


( 6 ) 


where c-y is the specific heat at constant volume and where C3 is a 
constant. The energy equation is 


uSx + vSy. = 0 


(7) 


For flows derived from Isentroplc Irrotational flows by a shock, there 
is the so-called condition of isoenergetic flow 




(7 - l)p 


= c. 


(5a) 


where 09 is a constant throu^out the flow, before and after the shocks. 
This equation is then Identical to eqi;iatlon (5); it implies eqmtlon (7) 
when equations (l) and (2) are satisfied. 

The botmdary condition at the solid body is the same as before; that 
is, the normal component of velocity must be zero. The condition at 
infinity, however, ’is no longer the same, because if the velocity were 
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constant at infinity the flow would again he Irrotational. The correct 
houndairy condition in this case is that the pressiore must he a constant 
at infinity. 


Since the flow in a part of the domain is rotational^ the velocity 
potential cp obviously cannot he- used and it is convenient to introduce 
in such cases the stream function i|f defined hy 



Equation (7) indicates that the entropy S is a constant along 
each streamline. The entropy S can therefore he written in the 
following form: 


S = cvf(t) (9) 

and the equation of state (6) as 

p = c^e^^^^pT (lO) 


With some calculation, the rotation of the flow may he written in terms 
of ■'Ir as follows: 


CO = Vjj - Uy 

= C3p7e^(^)f(t)/(7 - 1) 
VAEIATIOHAL PRINCIPLES 


( 11 ) 


Instead of studying the houndary-value problems as formulated in 
the preceding paragraphs, it is sometimes more convenient to study the 
associated variational problems, especially when the exact solutions of 
the differential equations are very difficult to obtain. In such cases. 
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approximate methods of solution such as the Rayleigh-Ritz method and 
the Galerkin method will he extremely useful. 

For irrotational flows, two variational principles were given by 
Bateman (reference 8), one in terms of the velocity potential and the 
other in terms of the stream fmction. The first one is more suitable 
for studying the flow passing arbitrary bodies and has been used by the 
senior author and his associates in computing such flows in refer- 
ences 1 to 7* In the present case, the second principle which is 
expressed in terms of I'' should be used. 

If the pressure is 


p = F(p) - pF'(p) (12) 

where F(p) is some function of p and the prime indicates differentia- 
tion, the variational integral is 


Jl = 


F(p) + 


(I'x^ *v^) 


2p 


dx dy 


(13) 


in which p and are to be varied Independently. 


Equation (12) is a differential equation of the Clalraut type. 
If p = c-]^p7, the solution of equation (12) is F(p) = C 2 P - 
The integral Jp then becomes 


7-1 


= 


CgP - 




7-1 


2p 


dx dy 


The condition that 6Jp = 0 then leads to 


^^2 ■ 


^1^ 

7-1 


p7-l _ _L, ^ 2 ^ 


2p 


2l-x 




+ (♦y/P)^S+ ax ay +/i g s* 


6p dx dy - 


ds = 0 


(11^) 
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where the last integral is taken on the houndaiy of the domain. If this 
line integral is zero, the condition = 0 leads to the following 

equations 


i?(’ 






C2 


ft/plx + (♦y/p)y ' 0 


as the Euler equations in the calculus of variation. The first equation 
is the Bernoulli equation and the second is the condition of irrotation- 
ality,. The condition of continuity is satisfied automatically by intro- 
ducing the stream function If. Th\is if the line integral is equal to 
zero the vanishing of the first variation of J leads to the desired 
differential equations. If the line integral is not zero, as in similar 
cases discussed in references 1, 2, 6, and 10, the value of the integral 
m\ast he subtracted from and sljould he modified accordingly. 

It is interesting to note that althou^ the variational integral (13) 
is formulated hy Bateman for irrotational flow only, it is also valid 
for rotational flow if in equation (12) F(p) is replaced by F(p,S) 
and is solved from the correct equation of state, equation (6) or (lO). 

In this case 


F - p ^ = C3ef(^)p7 

8p 


The solution of this differential equation may be written as 

c.e^(^)p’' 

F(p,S) = CpP ^ 

^ (7-1) 


(15) 


and a variational integral for the rotational flow after shock waves may 
be written as follows 


f 

/r< 



jj 1 

if . (7 - 1)J 

J 


dx dy 


(16) 
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This is the same integral recently obtained by Lin and Rublnov (refer- 
ence 9) • As before, the condition SJ 2 = 0 leads to 




ds*= 0 


The last Integral is again teiken on the boundary of the domain. If this 
integral is zero, there result 

* '*/) , rp. 

2p2 (7 - 1)P ^ 




as the Euler’s equations. The first equation is the energy equation (5a) 
and the second eqxiation is the rotation equation (ll). With the aid of 
the second equation, the equations of motion (l) can be derived from the 
first equation as follows: Writing the first equation in a sll^tly 

different form. 


— ^ + ^(u^ + v^) = C2 

Differentiating it with respect to x, 

7C2e^'(''^)p7-2p^ + - 7 - c^e^^^^pZ-lf ’(ilr)i^-^ + 


+ ^x = 0 


( 18 ) 
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Since p = 

Substitution in equation (l8) results in 


^ ^ ^ ^ *(t)i ^x ■'■ ■*■ ''^x ^ ° 

As C2P^e^(^)f‘(ilr)/(7 - l) = v^ - Uy and ^ tx = v, equation (I 9 ) is 

identical to the equation of motion in the x-direction. Similarly, 
differentiation of equation (l?) with respect to y resiilts in the other 
equation of motion. 

SHOCK COrJDITIOWS AS A RESULT OF VARIATIONAL PRINCIPLE 

In discussing the mathematical formulation of the problem, the flow 
as shown in figure 1 may be put into a simplified form as shown in fig- 
ure 2, where and Dg denote the regions when the flow is irrota- 

tional and rotational, respectively, and ^ 2.-2 shock wave plus 

the common boundary. The variational principle is that 

5J = 5J|j^ + 5 J 2 = 0 

leads to the desired differential equations in the corresponding domain. 
Here 
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where as shown In the above expressions is integrated over 

and Jg over D 2 ; Pi^''^! an<i P£^^2 that p and -are 

different functions in the domains and D 2 ; and the func- 
tions P 2 J and ^2 “to ^ varied independently. The 

condition 5J = 0 therefore leads to 


ff 

„ 7 CiPi^-l 


JjDi 

^ (7-1) 

2 Pl^ 


5p 


1 


dx dy - 


(W'"i)x ("‘"ly/'^i): 


dx dy + 




7C2e^(^2)p2^~^ 

(^2x^ + Tlf2y2j 

^2 

(7 - 1) 

2P2^ _ 


6 P 2 dx dy 


r 1 , 

n ^ ds + 

^^1 Pi 9n 1 


61^2 dx dy + 


Pi 
1 


1 ^1 


1 

'C, P, 3n ^ 1 5T "♦ 2 ' ^ " ° 


Consider the line integrals first. On a part of Cj and C 2 where 
the boundary is the svirface of the body, ilfj and ^2 equal to 

the chosen constants. Then = 51^2 " first two line 

integrals are zero. These integrals may not be zero when taken over 
the other parts of the boundaries. In such cases, the values of the 
integrals must be subtracted from 6 J and J must be modified accord- 
ingly. On the shock waves, if 




( 21 ) 
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and 


= ^2 ( 22 ) 


the third integral vanishes. Condition (21) requires the velocity 
components tangent to a shock wave to he the same in passing throu^ 
the shock. Condition (22) is the continuity equation in crossing a 
shock wave. On the common hotuadary, since it is a streamline^ 

= 6Tlr2 = 0 and the line integral vanishes. 

T^hen the line integrals are zero or modified to he zero, the 
following Euler equations are obtained from the well-kaown rules of the 
Calculus of Variation. They are 




(7 - l)pi 


2Pi 


(23) 


yP2 , 
{7 - l)p2 


2p, 


C2 ■ in Dg 




(24) 


(25) 


I'ac) - 1) In Da (26) 


Equations (23) and (24) indicate that the energy in the fluid remains 
unchanged in crossing a shock wave. With equations (23) eind (25) and 
equations (24) eind (26), the ecjuation of motion can he derived; that is, 
the momentm of the flow also remains the same in crossing a shock. The 
continuity equation is already given as eqmtion (22). Therefore, from 
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SUBSONIC AND TRANSONIC FLOWS PAST A CIRCULAR CYLINBER 
BEFORE OCCURIffiNCE OF SHOCK WAVES 


As a first step to tlie solution of transonic flows with shock waves, 
the flow Before the occurrence of shock waves will he considered. The 
variational integral J contains two variables, p and Numerical 

calculation may be carried out by following the Rayleigh-Ritz procedure. 
Unlike the cases considered in references 1 to 6 where the velocity 
potential q) is used, the labor Involved in the computation by taking If 
and p as variables in the Raylei^-Ritz method becomes excessively 
large. Instead, Galerkin’s method was found to be much more simple 
in this case. The application of Galerkin’s method to compressible- 
flow problems has been discussed in detail in references 3 and 4 and 
can be briefly oxrtllned as follows: First the variables ijr and p 

are written in the form of series which satisfy the boundary conditions 
but with vindetermined parameters, such as 






m n 


P ^o y 5 ■'^nn*^mn 

m n" 


(27) 


where A n,n and B^^ are the undetermined parameters. Then A^^ 
and Bmn are determined from the conditions 


t+x/p)x + (ty/p): 




dx dy = 0 


( 28 ) 





5p 

^®mn 


dx dy = 0 


(29)' 


With A ^y^ and Bp^y, determined, equations (27) give an approximate solu- 
tion to the problem. 


It was found, however, that near the limiting Mach number where 
shock waves are about to occur, the p-series converges rather slowly and 
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many parameters have to he nsed. The amomt of work involved in the 
integration and in the solution of the resulting simultaneous equations 
with many parameters again becomes very great. A modified Galerkin 
method has been tried and it was found possible to reduce much of the 
numerical work. The method is essentially as follows; Instead of 
determining Bmn by Galerkin ’s method, these parameters can be solved 

in terms of A mn from the Bernoulli equation: 



by the method of equal coefficients. The parameters A ^^ then become 
the only unknowns and they are detennined by the Galerkin method. 

In the case of two-dimensional flow past a circular cylinder- with 
unit radius the boundary conditions are at r = <»: 


P = Po 


and at r = 1 


^0 = 0 


Consider the case where the circxilatlon is zero. The flow must then 
be symmetrical with respect to both the x- and y-axls. The boundary 
conditions and the conditions of symmetry are satisfied if ''I'’ and p 
are assumed in series of the following forms: 


t = PqU 


(r - i) si. 6 . E g si. ne 


(31) 


and 



Z 


cos 2n0 

r2m 


P 


m=l n=0 


(32) 
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Taking two parameters in equation (3I) tiecomes 



Substituting equation ( 33 ) into equation ( 30 ) and equating the coeffi- 
cients of and cos nO, there result 
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^24 = i =12^ 


B = 0 if n ^ 5 
2n 


B 


3° ■ Mo2 


L_ 2Mo2(aii - Aii^ - 3 Ai3^) - (biiB2i + 


®31 m^2 


+ 3 A ^3 - 6 AiiAi 3 ^ - (2B12B2O + ®H ®22 + 

B12B2 I + Bi2B23)s + g " ~(i + i ®11®12^ 


= 


3= V - 1 


(®11®21 ^ ® 11®23 •" ^12®20 ■" ®12®2l()= ’■ — S ril “12 




h 3 

2 12 , 


B 


33 Mq2 . 1 


- 3 «oV - ( V 22 ^12^21}^ * * 


3®llV) 


B_, = 


3“ M„2-1 




=■ - 


35 -Mo^ - 1 




8 


11 12 


B.^ = 


36 M 2 _ 


Mo- - 1 


-®12®2lt® + Bj2 
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B3^ = 0 if n > 7 




1 


Mq2 _ 1 

2 


-f-( 5 Aii 2 + 9A132) - h(b2o^ + + JB232 + 

2 ^ 2k ^ + + B12B32) + + iBn%2 

JB 11 B 12 B 21 + Jb^^B^ 23 + i®12%0 + 5®12%4) + 

|®11^®12^ + i®12^) 


hi ~ .T2 


Mq- - 1 

B, 

Br 2 


¥(- 


+ 12A. 


11^13) " ®(‘ 


2B20B21 + ^21^22 ■*■ ®22®23 


B 


'42 


23®24 + 2B^iB30 + Bi1®32 + ®12B31 + ^12^33) + ^^^ (^^12^^21 

t®ll^®23 ®11®12®20 + ®11®12®23 I®11®12®24 2®12%l + 

-3Mo\iAy - ^ ^ ^ ^ ^ 


Mq2 - 1 


B 


11=33 * ^..=30 - =12=34) ; - ki%2 ^ Ki %4 - 

=11=12=21 |=iiBi^B23 + ^12%^) + . (7 ^ 1)7(7 - ^ 

?=11^=12^) 
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l[f (™20®23 + ®21»22 * \ t ? 2 \ * * 

2 B^B3 i * 2 B^B33) ^ - Ki %3 ^ 

®11®12®20 + ^11®12®22 + |BuBi2S24 + 

- ^> (|bu 3 b^ . 


m:2\ 


*‘‘5 ' 1^,2 

(7 


V - T[:“Vi*'22 ■" ^0®24 ®2*3 ^ ^11=33 + B11B35 + B12B33 + 

®12=3 s) ^ 2 ^'’'(i®l/%2 + ^ll\k * 1611^12^21 * ^1^12623 

^ 12^0 - |b^^B 2 ^) * l)j;(r - ^ ^ 

- i£®(^2X®21|. ®22®23 * ®11®34 * ®1I®36 * ^12^33) * 

i^{hl %3 * ¥11^1^22 * ^1^12*24 + ^12"62i) ^ 

- ^> (1^11^12 - ih ,^ J ) 

" =22624 - =11635 " =12=34) * H ^ ihi\k * 
^U®12=23 * ^12\3) * - * ' ^X i=ll^=12^) 

* ®12®35) * ■^'^1^(1=11=12=24 * ^12%3) ^ 


Iz. + 1 ) 7(7 - 2 ) ( 2 ^ 3 3 "\ 

2h- \^ 11-®12 ) 



KACA TN 2539 


19 


® 12 ® 36 ) 

(7 + 1)7(7 - 2 )/l^ 1 A 1 

^ If 12; 



+ 


(7 + 1)7 
2 



+ 


(31^) 


= 0 if n ^ 9 


y /y j_ T ^ P 

where H = — ' 7 For given free-stream Mach number and 7 , 

2 2(7 - 1) ■’ 

the *a are now expressed In terms of and ^ 

numerical example, 12 parameters are taken in the p-series which seems 
to give satisfactory convergence at high Mach numbers. 

The condition of ixrotationality is now the only differential 
equation remaining to be satisfied. In terms of \lf and polar coordinates, 
equation’ (h) becomes 

Pr'lt ° (35) 


where 


At = ^ ^ 


The Galerkln equation for determining the parameters A ^ti Is: 



+ 

r^ 


- p(At) 


Bt 


r d0 dr = 0 


( 36 ) 
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With the two parameters App and A^^^ ''1''^ eqiiatlon (36) becomes: 


pr=oop 0 = 2 jt 

/ / 

d r=i'^ e=o L 


p^Tjr + - p(At) 


sin 0 d0 dr = 0 (37) 


n r=oon e=2n 
^ r=l^ 0=0 


^r^r 


sia 36 a© dr = 0 (38) 


Carrying out the integration and taking 7 = 1.405 the following 
equations are obtained for different values of Mq. The parameters 

talffin in the p-series are 2 -]^-^, B-|p, ® 21 ^ ® 22 ^ ® 23 ^ ® 30 ^ 

B 31 , ®40-’ ® 4 l-’ ®42* these B^^ *s there result 

at Mq = 0,3: 

t 

0.7914i66Aii + 0.0148886A3_3 = -0. 0529632 + 0.0011372A^3A^3 - 

0.02903Q^2.2^ - 0 . 0353153 Aj_ 3 ^ + 

0.028i407Aii^a^3 - 0.0368037A2^^a^2^ " 

0 . 026577 Aii 3 _ o. 0537032 A 3^33 

-0.04i60-32A3^3^ + 2.1i13662A3^3 = 0.0565330 - 0.0524788a^^A^2 - 

o.oo 49465 Aj^^^ - 0 . 05 i 6919 A^ 2 ^ - 

o.o426662A^3^%^2 + 0. o6i4285A^^Aj^2^ + 

0.0182i44A3^3^3 _ 0. 0776527A^^3 
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At Mq = 0.1^: 
0. 8939718A3_3^ - 


-0. 1165173A11 


At Mq = 0,48: 
1.0588773A^^ - 


-0.2295625%! + 


0,029304-5A!2 = -0.1386222 - 0.0363370A A 

n 13 

0. 0208205A!!^ _ 0. 1594l51A.!2^ + 

0.0557317A!!% - 0.1400827A^A,_2 _ 

X-L Xj 

0. 0837686A!!^ - 0,i 685111A!2^ 
f 2.2168408 Ai3 = 0.1313933 - 0. 1462311Ai!A!3 - 

o.io92484A!!^ - o,2046769A!3^ - 

0.053295Aii2a^3 + 0. 3247757A!!A^3^ ^ 

0.0591811A!!3 _ o.1892021A!33 

0.0792806a^^ = -0.3375697 - 0.0768956a^!A^3 - 
o.03o879A^^2 - 0.7005271A^^^ + 

0. 0580232A!!2^^^ + 0. 362238A^^^A^^2 _ 

0. 0786318A!!3 _ 0. 19o4891A!33 

2.215112A!3 = 0.2795736 - 0. 7384032A^^A^^ - 
0.227209M!i2 _ 0. 9I22836A132 + 

0.210928A!!2^!3 + 1. 067842A!!A!32 + 
0.1269145A!i3 _ 0 . i446766A!33 
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At Mq = 0.5: 

1.0588045A11 - 0.09l^3964■A]^3 = -0.435^569 - 0. 07¥t03fiAiiAi3 - 

0.0687631Aii^ - 1 . 0282597Ai3^ - 

0.0279679Aii^Ai3 + 0. 390^25A3^3A^2^ 

1.236569Aii3 _ o.5792328Ai33 

-O.2586206A11 + 2. 151i-37^9Ai3 = 0.3i<-713l3 - 0.9668059A3^^A^3 - 

0.289680ifAii^ - 1.2930194 Ai3^ + 

0.4529877A3^]^2a3^3 + 1 . 7020295A3^3^A^2^ + 

0.136i6i7Aii3 „ 0. 3869221A3_2^ 


These equations can he solved hy the method of successive approxi- 
mation as outlined in references 1 and 2. The parameters at various 
Mach nTimhers and the computed maximum velocities at these Mach numbers 
are shown in table I. These values shew good agreement with those 
obtained by other approximate methods. At M^ = 0. 5 the method of 
successive approximation fails, and it was foimd by graphical method 
that these equations do not have a solution. This agrees with a previous 
investigation as reported in reference 7 where the velocity potential cp 
was used. This indicates that the flow is no longer Irrotatlonal and 
shock waves have probably appeeired. The limiting free-streeim Mach 
number is between O.W and 0.5. This value is probably too hi^ because 
only two parameters in t are used. This fact has been discussed in 
reference 7- 

To get an idea of the convergence of the and p-series, the values 
of maxirnm velocities and p at Mq = 0.48 are computed and listed in 

tables II, III, and IV. ' The p-series does not appear to be convergent 
when the parameters are added one by one, but if these parameters 

are added by the group the series becomes convergent. This can be seen 
by comparing the values of p and the max imum velocities with various 
numbers of Bum’s as given in tables II and III. 
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TBAKSONIC FLOW PAST A CIRCULAR CYLDTOER WITH SHOCK WAVES 


When the free-stream Mach mamher reaches 0. 5 there is no solution 
to the equations obtained by assuming irrotational flow in the entire 
domain. Shock waves probably have occurred. Thus the flow is Irrota- 
tional only in the region before the shock waves or and becomes 

rotational in the region behind the shock waves D2. The approximate 
solution according to the modified Galerkin method can be carried out 
as follows: 

(1) The locations and the obliqtiity of the shock waves are first 
assumed from the best evidences available. 

(2) When the assumed shock waves are Inserted in the flow, the 

domains Dp and D2 are thus defined. In the region the flow 

is still irrotational. The boundary conditions are that the velocity 
at infinity is eqtial to the undisturbed velocity and on the surface of 
the body must be equal to a chosen constant. The values of 

and may be assumed in the same form as in the case before the shock 

waves occur. The parameters Ajjjjj and are determined from the 

modified Galerkin method. The only difference between the present case 
and the previous case is that in this case the integrals are to be 
extended only in the region instead of the entire flow region. 

(3) The entropy distribution S = c^f(v) in the flow after the 
shock is to be assumed next. From the flow conditions before the occur- 
rence of shock waves, a good approximation of f('»J'’) may be found by 
rou^ preliminary calculations. 

( 1 )-) The boundary condition at infinity of the flow after the shock 
is that p is equal to the constant pressure Pq. With f(''k) assumed, 

Po(x^y) can be computed from equation (lO) and the velocity distribution 
from equation ( 5 a). 

(5) With the above boundary conditions and the condition that 

on the assmed shock waves and common streamlines" of and I>2y "^2 

and P2 are then assumed in the forms of series with undetermined param- 
eters. These parameters can again be determined by the modified Galer kin 
method. 

( 6 ) After these undetermined constants in and ^2 computed, 

the values of (^/pl) (^^l/^°) (^^2/^) shock waves 
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can "be found and 
to ^l/P2^ 


is calculated. If 
f('*li’) foimd is not 


(Vpi) (^l/^'^ equal 

the same as fCI^) assumed. 


as is generally the case, adjust the locations and the obliquity of the 
shock waves. Then based on the new position and shape of the shock waves 
assimie a new f(‘'l'^) and repeat the computations again until these values 
check with each other in two consecutive cycles. 


} 


A numerical example will now be carried out. With and as 

assumed ih equations (31) and (32) and by taking 2 parameters App 
and Ap3 in equation (31) and 12 parameters in equation (32), it was 
found that there does not exist a solution at Mq = 0.5. In refer- 
ence II it is shown that when the shock waves first occur they are 
almost normal shocks located near the point of maximum velocity. As 
a first approximation, assume that at Mq = 0. 5 there are two normal 
shock waves located at 0 = it/2 and 3Jt/2 from r = 1 to r = 1.5. 

The separating streamlines between the regions Dp and I>2 are assimied 
to be the same as the streamlines -in the flow at Mq = 0.14-8 passing 
through the points r = 1.5 and 6 = it/2, 3ji/2. 

To determine the parameters Ap^^ and A-j^^, integrals (37) and (38) 

eire to be extended in the above assimied region. This can be done by 
taking first the integration over the entire flow region and then sub- 
tracting from the resvilting values the valixes of these integrals evalu- 
ated in the region D2. The region D2 consists of two areas, namely, 

a semicircular ring with r from 1 to 1.5 and 0 from ir/2 to -rt/2 and 
a symmetrical tail-shaped area. The dominating parts of the values of 
the integrals happen to be in the semicircular-ring area where the 
integration can be carried out analytically. In the tail-shaped area, 
analytical integration becomes rather laborious and numerical integration 
has been ■used. After subtraction, these two equations are 

2.0587564A11 - 0.2021789A13 = -0.8773276 - 0. 1912757AX1A13 - 


0. 1931520Axx^ - 2.010768Ax3^ + 

0.3ij-90648A3^X^%3 + 0.891^8537A^^A^2^ - 


0.0087952AX1 + 9.21901h8Ax33 
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-0.335855^Ai]_ + 5.296336IA13 = 0.76527^8 - 1. 747661 <- 7 AiiA 3^3 - 

0.4624384A;j^j_2 _ 2.088238^13^ + 
0.7361793Aii2a^2 + 2.903363U11A132 + 
0.3574438Aii3 - 7.1979^^3^133 

While a solution does not exist before the integrals in D2 sub- 

tracted, solution of these equations now does exist. Solving, there 
result: 


All = -0.437 A13 = 0.100 

With these values of An 'and A13, the point at 6 = nf 2 with a 

unit local Mach number is computed to be at r = 1. 6. This indicates 
that the shock waves probably would terminate at a value of r between 
1.5 and 1.6, and the region Di as assumed is a good approximation. 

The boundary conditions for ilf2 and P2 are as follows: 

At r = to, p = Po and therefore 

02 = if (39) 

At r = 1, since iITi is taken as zero, 

^2 = 0 (^0) 

To insTire the vanishing of normal velocity, 

'1^0 = 0 (4l) 

On the shock waves and the common boundary. 


m 


■^2 = ^1 
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Conditions (39) j (^0)^ and (4l) are satisfied if -^2 Pg 

assumed in the following forms: 



and 


P2 - 



-f 

e 


{^2) + 

r=oo 


B', 


m=l n=0 


mn 


cos 2n0 


(H) 


Note that ''^2 becomes a function of y = r sin Q at r = <». The 
constants are to he determined hy the collation method in terms 

of so that ~ iiuniBer of points on the shock waves 

and the common boundary. With ^('* 1 ^ 2 ) assumed, B'nm solved 

in terms of A’^^ from equation (2h) and the parameters A’jjjj^ may 

he calculated from equation (26) hy using Galerkin’s method. A numeri- 
cal example has been tried. It was found, however, that imless some 
modern hi^-speed computing machine is used the computation becomes 
very lengthy. 


New York University 
New York, N. Y. , 


March 1, 1951 
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TABLE I 

PARAMETERS AM) COMPUTED MAXIMUM VELOCITIES 
AT VARIOUS MACH NUMBERS 


■" 1 

0.3 

o.k 

0 . 11-7 

0.48 

^11 

-O.O6TT 

-0.1537 

-0. 296 

-0.316 

A13 

0,0255 

0 . 0 ll -04 

0.088 

0.087 

p . ^ 

5 - 

0 = 90° 

0 . 825 ^^ 

0. 6972 

0.3705 

0.2680 

<1 

TT 

^ 0 = 90° 

2.1972 

2.2859 

3. ^870 

4.5525 
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TABIE II 


CONVERGENCE OF p-SERIES AT = O.i^ 

1 


Parameters 
taken in 
p series 



"13 


r = 1 
0 = 90 ° 


r = 1 


Bll 

®11®12 

®11®12 

®20 

®11®12 

®20®21 

®11B12 

B20B2I 

B22B23 

%1®12 

®20®21 

B22B23 

B30B31 

B32 

®11®12 

B20B2I 

B22B23 

B30B31 

B32B40 

B4 iB 42 

-0.0988 

-0.0836 

-0.490 

-0.390 

-0.3950 

-0. 3547 

-0. 316 

0.0777 

0.0781|- 

0.120 

O.O9I' 

0. 0960 

0. 1110 

0.074 

0.69^3 

0. 6522 

0. 14l6 

0. 2931 

0.2686 

0.2788 

0.2680 

2.3765 

2.5698 

5.5104 

3.5416 

3.7902 

4.5503 

4.5525 
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^ TABLE III 

CfflJVERGSHCE OF p -SER I ES AD? M 


Humber of 
paraiaeters 
taken In 
p-serleB 
( 1 ) 


r = 1 


^0 9 ^ 90 ° 


r = 1 


0 = 90° 


o 7Yf;'^l2-'i6Q8l‘5.510U 


5 

6 

d 

0.2558 

3.7902 

3.9797 


0. 691^3 1 0.6522 1 0.l4l6l 0.2931 0. 26861 0. 2558 0. 5187 1 0.38531 0.27881 0, 0378 lo.i^520l 0.2680 


11 

12 

0.4520 

0.2680 

2.6966 

4.5525 


^Parameters In p-aerles are taken according to the following order: 

®23^ ®30^ ^1» ^32^ ®4l^ ®42* 

































NACA TN 2539 


31 


TABLE IV 

CONVEEGMCE OF t SERIES WITH 
12 PARAMETERS IN p-SERIES AT Mq = 0.46 


Parameters taken 
in ■'If-series 

All 

■S.1^13 

r = 1 

^ at 

8 = 90° 

0. 1007 

0.2680 

q "■ = 1 

— at 

U 0 = 90° 

13.5194 

4.5525 
















